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We employ an approach wherein vacuum entanglement is directly probed in a controlled manner. 
The approach consists of having a pair of initially nonentangled detectors locally interact with 
the field for a finite duration, such that the two detectors remain causally disconnected, and then 
analyzing the resulting detector mixed state. It is demonstrated that the correlations between 
arbitrarily far-apart regions of the vacuum of a relativistic free scalar field cannot be reproduced 
by a local hidden-variable model, and that as a function of the distance L between the regions, the 
entanglement decreases at a slower rate than ~ exp( — (L/cT) 3 ). 



The vacuum state of a relativistic free field is entan- 
gled. For two complementary regions of spacetime, such 
as x < and x > 0, this entanglement is closely re- 
lated to the Unruh acceleration radiation effect ID, and 
gives rise to a violation of Bell's inequalities 0, S- For 
two fully separated regions, entanglement persists |j|, al- 
though, it is not known whether Bell's inequalities are 
violated, or how entanglement decays with the increase 
of separation, as compared to correlations. Similar ques- 
tions concerning entanglement_have been addressed in 
the case of discrete models 
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In this letter we shall study this problem by probing 
the field's entanglement with a pair of localized two-level 
detectors [Sj. This is done as follows. A state is prepared 
in which the two detectors are not entangled with one 
another, or the field. We then have each of the detec- 
tors is made to locally interact with the field for a finite 
duration, such that the detectors remain causally discon- 
nected throughout the process (Fig. Since entangle- 
ment cannot be produced locally |2j , once the interaction 
is over, the net entanglement between the detectors must 
necessarily have its origin in vacuum correlations. The 
interaction thus serves as a means of redistributing en- 
tanglement between the field and the detectors. We shall 
show that for arbitrarily far-apart regions, the detectors' 
final mixed state, after filtering, violates Bell's inequal- 
ities, and in the process obtain a lower bound on the 
amount of vacuum entanglement. 

To set-up the model, we shall assume that the detec- 
tors are localized within a region of a typical scale of 
i?, and are separated by a much larger distance L » R. 
Consistency with relativity requires us to use detectors of 
a rest-mass M, for which R» Xcompton = jj^- In this 
limit, the effects of both detector pair-creation, and the 
"leakage" of each detector's wavefunction to the outside 
of its localization region, become exponentially small, of 
the order of - exp(-2^) [13, E|- Note, that this en- 
sures that the overlap between the detectors' wavefunc- 
tions is negligible. Under these conditions, in their rest 
frame, the detectors can be described as nonrelativistic 
quantum mechanical systems. Finally, we shall assume 



FIG. 1: The world lines of detectors A and B are shown for 
the duration of the interaction. The horizontal and vertical 
axes are space and time respectively. The arrows denote the 
emitted radiation. Notice that the radiation emitted by de- 
tector A (B) does not affect detector B (A), since for t > T 
the interaction is switched-off. 



that, by means of an external coupler, each detector's 
degrees of freedom can be coupled "at will" to the field. 
Since the coupler need not be of the same type as the 
studied field, we shall make the additional assumption 
that it can be described classically, and therefore does 
not generate entanglement. 

There have been several proposals for detector models 
which can satisfy the above requirements; notably, the 
Unruh- Wald "particle in a box" detector 0, and the 
DeWitt monopole detector model 0|. In both models 
the detector Hamiltonian is §0"z, with being the energy 
gap between the two levels and a z a Pauli matrix. The 
field-detector interaction Hamiltonian is 
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4>(x, t) is a relativistic free scalar field in three spatial 
dimensions, the a are the detector's energy raising and 
lowering operators, and e(t) governs the strength and du- 
ration of the interaction. The function ip(x) is a function 
of the detector's spatial degrees of freedom, and is deter- 
mined by the model employed 0, ^| . 
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FIG. 2: Emission and Exchange Processes. The exchange 
amplitude, portrayed on the right, is dominated by a single 
off-shell emission followed by an on-shell absorption, while 
the double emission amplitude, portrayed on the left, consists 
of two off-resonance processes. Thus, as Qa,b increases, the 
emission term, \\E A \\ \\Eb ||, decreases more rapidly than the 
exchange term, |(0|Xas)|. 



\ J J dtdt'f[H A (t)H A (t')] + {A^ B). We observe that 
in the first term above, the state of the detectors is un- 
changed, while in the second term both detectors are ex- 
cited and the final state of the field is \Xab) = ^Ji'&ijIO). 
Since \Xab) contains either two photons or none, it de- 
scribes, respectively, an emission of a photon by each of 
the detectors or an exchange of a single virtual-photon 
between them (Fig. |3J). Finally, the last couple of terms 
describe an emission of a single photon by either detec- 
tor A oy B. In this case the final state of the field is 
\E A ) = *+|0>, or \E B ) = *+|0). 

Tracing over the field degrees of freedom, when work- 
ing in the basis {|i)}| =0 = {TT, U, IT, U} and employing 
the notation H-XabH 2 = (Xar\Xar), we obtain the de- 
tectors' reduced density matrix [lfj 



Consider now a pair of DeWitt monopole detectors, 
A and B, that are localized about the coordinates x A 
and Xb, respectively. These detectors interact with the 
field through H int = H A + Hb, where H A and Hb are 
interaction Hamiltonians of the form of Eq. Q). The 
window functions e A (t) and esit) are chosen to vanish 
except for a finite duration T, such that cT << L — 
\xb — x A \, ensuring that the detectors remain causally 
disconnected throughout the interaction. In the following 
we shall work in the Dirac interaction representation and 
employ "natural" units (h = c= 1). 

Since the interaction takes place in two causally discon- 
nected regions, the Hamiltonians H A and Hb commute. 
The evolution operator U for the whole system thus fac- 
tors to a product of local unitary transformations 



U = Tie 



-i f H A (t)dt 



-i J H B (t')dt'-i 
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where T denotes time ordering. This guarantees that 
U does not change the net entanglement between the 
regions. 

We take the initial state of the detectors and the field 
to be I*.,) = | |a) | !b)|0), where | J.) and |0) denote 
detector and field ground states respectively. In the weak 
coupling limit ti{t) « 1 (i = expanding to the 

second order we get 



I*/) = 



(i-C)| U)-$+ $+| TT) 



- t$+l B \n)-rtA<S>B~\m\\0)+O(e 3 ), (3) 
where $f = / dte i (t)e ±<n * t / d 3 a;^(f)0(x, t), and C = 
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Note the two types of off-diagonal terms. The am- 
plitude (0\X A b) acts to maintain coherence between 
I IaIb) and | TaTb)j while the amplitude (E a \Eb) acts 
to maintain coherence between | [ a ]b) and | ]a[b)- It 
is the relative magnitude of these off-diagonal terms, as 
compared to the diagonal decoherence terms, that deter- 
mines whether the density matrix is entangled. 

A density matrix is said to be inseparable or entangled 
iff it cannot be expressed as a convex sum of local density 
matrices ll 71 . In the present case of a 2 x 2 system, a nec- 
essary |l 8| and sufficient 0] condition for inseparability 
is that the negativity [13] be positive. We shall therefore 
use the negativity as a measure of entanglement. The 
following expression is obtained for the negativity 



M( P )^mX A B)\-\\E A \\\\E B \\ >0. 
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Physically speaking, the inequality above is satisfied if 
the single virtual-photon exchange process is more prob- 
able than the off-resonance emissions of a single pho- 
ton by each of the detectors. The main contribution 
to the entanglement then arises from states of the form 
a\ UIb)+0\ UTb). 

The inequality (Eq. [SJ can be reexpressed as 
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where the factor e~ w R accounts for the "smearing" of 
the detectors, and €i(u>) denotes the Fourier transform of 
e.;(t). (Note that for a massive field a factor of , == . 
must be added to each integral.) The term sin(wL) 
on the left-hand side can be interpreted as an effective 
window-function that characterizes the field's response 
to the detectors. When integrating over ui, this func- 
tion governs the overall sign of each mode's contribution, 
and thus acts to reduce the exchange amplitude. This de- 
structive interference effect can be minimized by employ- 
ing a window-function ca for which sin(wL)e j 4(f7 J 4 + lo) 
remains positive over a finite integration regime in the 
limit of large L. A superoscillating function meets this 
requirement [2lt Eif . In particular the function [2^ 



e A (u) = /MJ VW-JV 2 ((W-1) 



(7) 



where f(to) is any function that converges faster than l/u> 
and has finite temporal support. We observe that 
is bounded in time as required, and oscillates like sin(wL) 
about u) = uj s ± >/~N /2L, where Tuj s — NL/T, approx- 
imately y/~N times, before gradually resuming "normal" 
slow oscillations for larger values of u. The use of su- 
peroscillations, however, is not without a price. For 
loT < N^J (L/T) 2 — 1 the function I a (w) decays expo- 
nentially, rendering the exchange term, and hence the 
negativity, exponentially small in L. The second window- 
function, £b, is a fixed hat function, convolved k times 

with itself. In uj space it assumes the form ( S%n ^Jl k " > J ■ 
For given values of L and T we choose the energy gaps 
n A = yVWtT ~ 1 and Mb = los - tt A » N/2L. This 
choice of Qb fixes the center of the window-function, 
£b{Qb — w), in the region of superoscillations. For large 
values of L, the inequality (Eq. EJ can then be approxi- 

> 1. This ratio can 



mated by V^(B 2 TOTt(^. 
be made arbitrarily large, at the expense of reducing the 

negativity, by increasing N. To get a lower bound, we 

take N = (^) 2 . Substituting the expressions for Qa and 

£Ib into the above approximation, it takes on the form 

(T/Lf(L/2T) k > 1. For k > 5 this ratio increases with 

L. We thus get a lower bound on the negativity 



-{L/cTf 



(8) 



Numerical computations show that this bound can be 
further improved. Taking N > L/T, we get M(p) > 
e -( L / T ) |24j|. The leakage of each detector's wavefunc- 
tion to the outside of their localization regions introduces 
a correction of the order of e~ 2MR to the above expres- 
sion. However, this correction can be made arbitrarily 
small by setting e ~^ L / T) ' » er 2MR . Note that we are 
free to do this, since the mass scale, M, and the distance 
scale, L, are independent. 

In passing, we would like to point out that in the case of 
the electromagnetic field or any other free field, the anal- 
ysis can easily be repeated. Similar results are obtained 




FIG. 3: Violation of the CHSH inequality. The dashed 
line represents the negativity of the detectors before passing 
through the filter, while the solid line represents the quantity 
M(p) — 1, which is calculated after the detectors have passed 
through the filter. Since the CHSH inequality can be written 
in the form M(p) < 1 we see that for fl -> 00, the CHSH 
inequality is maximally violated. 



in the case of the finite duration coupling of the detec- 
tor's magnetic moment, or electric dipole to the field. For 
a massive field, in the limit of large L, the above result 
remains unchanged, because then the contribution to the 
integrals arises from the range of frequencies u> ^> m for 
which the field effectively behaves as if massless. 

The reduced density matrix derived in the previous 
section is entangled. The question arises as to whether 
these vacuum correlations admit a local hidden-variable 
(LHV) description [3} . Applying the Horodecki theorem 
[25j to Eq. 0J, we find that the detectors' final state 
does not violate the CHSH inequality 1261 . We shall now 
demonstrate that by using local filters J27j , a violation of 
the CHSH inequality can be achieved for every separation 
distance, L. Hence, the density matrix (Eq. 0} reveals a 
"hidden" nonlocality [3] , as in the case of Werner states 

H 

To show this we follow Gisin [2j| ■ Once the interaction 
with the field has been switched-off we have each detector 
pass through the filter 



fA,B 



I 

rj 
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The density matrix is thus transformed according to p — > 
Pf = {f A <H) /b)p(/a <8 /b), so that the (nonnormalized) 
filtered density matrix is given by 



II^abII 2 

rf\\E A f 

n 2 {E A \E E 

\ - v 2 (x AB \o) 



-r f 2 <0|X AB ) ^ 

n 2 (E B \E A ) 
'<7 2 ||£ B || 2 

„ 4 (i-||i? B || 2 -||£: A || 2 ) / 

Consider now the choice i] 2 1=3 (Xab|0). We note 
that the 00, 33, 03, 30 terms (Eq. I10[) are now nearly 
equal, and of the order of KX^slO)! 2 . Ideally, in the 
absence of decoherence terms in the inner block, these 
terms come close to reproducing a maximally entangled 
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state | IT) — | II)- Notice, however, that the decoherence 
terms are of the order of |(^'ab|0)|||£' j 4,b|| 2 - Previously 
we have shown that the ratio |(Y^ B |0)|/ (H-EUll H-E-sll) can 
be made, by a suitable choice of window-functions, arbi- 
trarily large. Therefore in this extreme limit the relative 
strength of the decoherence terms, as compared to the en- 
tangling terms, is greatly reduced, and /?/ can be brought 
as close as we like to a pure, maximally entangled state. 
This implies a maximal violation of Bell's inequalities for 
the hnal state of the detectors, and since initially the 
detectors are not entangled (thus admitting a LHV de- 
scription), it follows that correlations between arbitrarily 
far-apart regions of the vacuum cannot be ascribed to a 
LHV model. 

Maximal violation can be achieved at the price of re- 
ducing the detectors' entanglement (negativity), which 
grows smaller in the above limit (Fig. |3J). We now wish 
to quantify the more general case, for which the final 
state is more entangled (larger negativity), but gives rise 
to a weaker violation of Bell's inequalities. Applying yet 
again the Horodecki theorem , we find that the CHSH 
inequality is violated iff 

\(0\X AB }\ \\X AB \\ 
\\E A \\\\E B \\ \(0\X AB )\- [ ' 

Note that for-apart regions, this inequality is only slightly 
stronger than the condition for entanglement. 

Interestingly, if we repeat our process many times, the 
resulting ensemble of entangled pairs can be reduced to 
a smaller one of higher quality entangled pairs. This 
process is known as distillation of entanglement, and is 
feasible for any inseparable 2x2 mixed state 30] . Fur- 
thermore, since Bell's inequalities are violated in our ex- 
ample, the two detectors could be used directly for tele- 
portation tasks, without having to distill them first |3lj . 

In conclusion, we have presented a new physical effect 
of vacuum fluctuations which is associated with quan- 
tum nonlocality. This effect stands in marked contrast 
to other vacuum phenomena, such as the Lamb shift or 
the Casimir effect, which to some extent can be "mim- 
icked" by classical stochastic local noise |32j . 
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